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Abstract 

We analyze patterns of remnant discrete symmetries that arise from \J(1) N 
theories by spontaneous breaking. We describe a simple, geometrical way to un- 
derstand these patterns and provide methods for identifying the discrete symme- 
tries and bringing them to the simplest possible form. Applications in GUT and 
string model building are briefly discussed. 
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1 Introduction 



Symmetries play a key role in our understanding of fundamental physics. While forces 
originate from continuous symmetries, discrete symmetries turn out to explain many im- 
portant properties of matter. The perhaps most fundamental examples are the reflection 
and conjugation symmetries P, T and C, yet there are other crucial 7L^ symmetries. 
Examples for such symmetries include matter or R parity in the minimal supersymmet- 
ric extension of the standard model (MSSM) [1]. Matter parity (or family reflection 
symmetry [2]) leads to a suppression of the proton decay rate and explains the stability 
of the MSSM dark matter candidate, the lightest supersymmetric particle (LSP). It can 
discrete 7L>2 subgroup of a baryon-minus-lepton-number symmetry U(1)b-l- 
One can break U(1)b-l to matter parity by giving vacuum expectation values (VEVs) 
to fields with B — L charge ±2. In SO(10) grand unified theories (GUTs) one may switch 
on appropriate components of a 126-plet (see e.g. [3]); in string theory an analogous 
breaking pattern can be achieved [4]. Discrete symmetries might also be an important 
ingredient for solving the flavor puzzle. 

Obtaining discrete symmetries as remnants of gauge symmetries is, at a fundamen- 
tal level, motivated by anomaly considerations. While discrete symmetries often are 
imposed ad hoc as global symmetries, it has been argued that global symmetries are vi- 
olated by quantum effects unless they originate from a gauge symmetry via spontaneous 
symmetry breaking [5-7]. 

Here we focus on Abelian discrete symmetries i.e. Jjn groups. (For a recent discussion 
of how to obtain non- Abelian discrete symmetries by spontaneous breaking see [8].) It 
is well known how to obtain a single Zjv-symmetry from a U(l); this will be reviewed 
below. On the other hand, the most promising candidate for a consistent theory of 
quantum gravity, string theory, typically provides us with models which exhibit a large 
rank gauge symmetry. This symmetry has to be (spontaneously) broken to the standard 
model in realistic vacua, which will generically lead to a non-trivial set of remnant 
discrete symmetries. The purpose of this study is to work out how such symmetries can 
be identified and put into a simple, i.e. canonical, form. 

2 Multiple %n symmetries from several U(l) factors 

In this section we show how to determine the remnant symmetries if we break multiple 
U(l) factors by a set of VEVs. Throughout this study we shall assume the U(l) charges 
of all fields to be integers, which implicitly fixes our conventions for charge normalization. 
We start with the easiest example of one U(l) and two fields. 



Consider a theory with gauge group U(l) and two complex scalar fields and ip (cf. [5]). 
Under the U(l) the fields transform according to 



2.1 Review: U(l) 7L q 
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i) -> e -' la{x) iJ (lb) 

with q G N, i.e. has charge q and ^ has charge —1. The terms 

0*0, , 0*0-0> and 0V 9 + h.c. (2) 

and powers as well as products thereof are gauge invariant. Suppose now that acquires 
a VEV. This leaves us with effective interaction terms of the form (■0 9 ) n with n G N, 
dictated by the symmetry 

V> -> e 27rU/9 V with £ = 0,1, ...g-1 . (3) 

An equivalent way to obtain this result is stating that the remaining symmetry is deter- 
mined by the condition 

e iqa ^(j) = =>. qa = 2n£ with feZ. (4a) 

Hence, 

|^ = modi, (4b) 

27T 

or, equivalently, a G — Z. We have just rederived the well known result that, by giving 
a VEV to a field with charge q, the U(l) gets broken to a Z g discrete subgroup. In 
what follows we will generalize this to situations in which several U(l)s get broken to a 
number of Z n s. 



2.2 The general case 

Let us now consider the general case of a \J(l) N gauge theory with M scalar fields 0^ 
(1 < i < M), which will acquire VEVs, and K other 'matter' fields ^ (1 < j < K). 
We will denote the charge of the fields w.r.t. the j th U(l) factor by qj{4>^) and qj(ip®) 
respectively. Accordingly, the 0^ fields transform as 

0« - exp^i^^W)^.^ 0W. (5) 

q(<f)^) can be thought of as an A/"- dimensional charge vector and = ^(0^) as an 
M x N charge matrix. Suppose now that iV > rankQ. In this case, there are unbroken 
U(l) factors. Then we can rotate the U(l) directions by an orthogonal transformation 
such that all 0W will be uncharged under (N — rankQ) U(l) factors. These U(l) factors 
will not be affected by the VEVs of the 0W fields and we do not have to consider them 
any further. Therefore, without loss of generality, we will from now on consider the case 
N < rankQ. Notice also that in supersymmetric theories the rank of the charge matrix 
cannot be maximal as D-flatness requires a non-trivial solution of s ^ Ji niq{4>^) = 
with rii G No- 
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To identify the remnant discrete symmetries after spontaneous symmetry breaking, 
consider the generalization of equation (jlj) in our simple example, 




i$>(* W )«i) ^ = <P (l) ■ (6a) 



This is equivalent to 

X)?»(^ (0 )«i = 27r£i with ^ez. (6b) 



As in (I4bl) . the right-hand side represents the usual 'mod conditions' for discrete break- 
ing. Equation (16bj) can be recast in matrix notation, 



Qa = 2n£ with £ e Z M . (6c) 

Recall that Q is an M x iV matrix with elements in Z. Such a matrix can always be 
brought into diagonal form by two unimodular transformations, i.e. invertible matrices 
over Z. Concretely, there exist A G GL(M, Z) and B G GL(JV, Z) such that 

AQi? = D = diag'(di, . . . , d N ) and di divides d i+ i . (7) 

Here diag' means that D is an M x iV matrix whose only non-zero elements di are at 
the ii positions. There are N non-zero diagonal entries due to rankQ > N. D is called 
Smith normal form (or just 'normal form' as in [9]) of Q. Note that unimodular matrices 
have two important properties: 

1. They have determinant ±10 

2. The greatest common divisor of all matrix elements in a single row is 1 because 
otherwise the determinant would not be ±1. The same applies to each column. 

Transformation ([7j) allows us to rewrite (J6"c|) . 

A' 1 DB' 1 a = 2ir£. (8) 

Now multiply this equation by A. Due to the second property of unimodular matrices, 
the 'mod conditions' in equation (l6"c|) remain unchanged, since £' = Ai still takes all 
values in Z A/ if £ does. Defining a' = B^ 1 a we arrive at 



p'j 

a'j = 2tt— with < £' J < d 5 - 1 . (9) 

dj 



x If the determinant was zero, the matrices would not be invertible. If the absolute value of the 
determinant was greater than 1, the inverse would not be an integer matrix, i.e. the matrix would not 
be invertible over TL. 
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Hence, we see that the remnant discrete symmetry is Z^ x . . . x 7L& N . If there are some 
di = 1, the corresponding factors are trivial and can be omitted. The fields ip^ then 
transform according to 

^ - ex P (iJ2Qk^ U) )B k£ a'A^ = exp ( 2vri ^ Qk(^ j) ) K~ J (10) 



Ar.f 



with new charges q'^ip^) = J^iQii^^) Afe; which are defined modulo c? fc . That is, we 
can choose q' k (ip ij) ) G {0, . . . , d k - 1}. 



2.3 An example with two U(l) factors 

Let us illustrate the above procedure by an example. Consider a U(l) x U(l)' theory 
with three fields obtaining VEVs and two other fields. That is, we have iV = 2, M = 3 
and K = 2. The charges are given in table [TJ In this example, we only consider scalar 



(a) VEV fields. 

U(l) U(l)' 
0« 8 -2 
0( 2 ) 4 2 
0( 3 ) 2 4 



(b) Matter fields. 

U(l) U(l)' 
1 3 
^ 1 5 



Table 1: Charges of the fields with respect to the two U(l) factors. 

fields, such that we do not have to worry about anomalies. Later, in the applications in 
section [5] we will discuss supersymmetric, anomaly-free settings. The charge matrix is 
given by the charges of the VEV fields (cf. tabled] (a)), 

/8 -2\ 

Q = I 4 2 I . (11) 

The diagonal matrix D and the transformation matrix B are 

D = I 6 1 and B = ( J " 2 ) . (12) 

Hence, we can read off that we are left with a ^2 x ^ symmetry. The charges of the 
fields can be determined by multiplying their charge matrix by B from the right, 




1 3 W 1 -2 \ = ( 1 1 
1 5 ) \ 1 I 1 3 



(13) 



The new charges of the ip 1 ^ fields are given by the rows of this matrix. Altogether we 
find that the setting discussed here leads to a Z2 x Z@ symmetry, under which ip^ 1 ' has 
charge (1,1) and ip^ has charge (1,3). 
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2.4 Visualization 



We will now provide a simple, geometrical way of envisaging the symmetry breaking 
patterns. First, notice that a direct product of groups x . . . x Th& N can be represented 
by an iV-dimensional lattice. Each element of the group can be thought of as a point 
in the fundamental region or unit cell of the lattice. The volume of the unit cell is the 
number of elements and the order of the group. 

Let us illustrate this by the above example. The VEV fields <^> with charges as listed 
in table [T] (a) span the lattice. Since the first VEV field is an integer linear combination 
of the latter two, = 3(p^ - 2<p^\ ^ and 0® span a basis of the charge lattice as 

The matter fields if)^ are represented by the 



illustrated by the arrows in figure 
bullets. A coupling (^ (2) )' 
if ni g(V 
squares. 



la 



is allowed by the discrete symmetries if and only 
+ n.2 <7(/0 ) hes on a node in the charge lattice, which are represented by 




V,(2) 



(b) Diagonalized lattice. 



(a) Original lattice. 



Figure 1: Illustration of discrete breaking. 



Our procedure described in section 12.31 amounts to finding an orthogonal basis for 
the lattice which is given by the rows of the diagonal matrix D. The matrix A performs 
a rotation of the ^ charges which eliminates linear dependencies, while B is the trans- 
formation between the bases. In the new basis, the lattice is orthogonal, and the ip( 1 ' 
charges are given by the projections on the basis vectors. We see that there is a Z2 x 7L§ 
discrete symmetry where the Z2 corresponds to the horizontal and the 7>q to the vertical 
direction in figure 1(b). The matter fields have charges (1,1) and (1,3), respectively, 
which are of course nothing but the coordinates of the bullets in figure 1(b) Note, 



however, that the true symmetry is smaller than that; we will describe in section |3] how 
to eliminate potential redundancies. 



2.5 Inverting the problem 

Often one is interested in 'inverting' the above procedure. Instead of determining a dis- 
crete symmetry that is left after certain fields attain VEVs, one would like to understand 



5 



A = 1 Vn a A (a) ; n a eZ} (14) 



which fields need to be switched on in order to obtain a desirable discrete symmetry. 
Examples for such symmetries include matter parity and proton hexality [10]. We de- 
scribe now a simple algorithm for accomplishing this task. We focus on the case of a 
remnant Z2 symmetry; the extension to the general case is straightforward. 

© Start with a set 'matter' fields tp^ with charge vectors q(ift( 1 '). 
® The q{ip^) span a lattice 

^n a A (a) ; n a G zj 

such that q{^ 1 ') G A, i.e. 

qty®) = ^n«A (a) with b}gZ 

a 

and there is no coarser lattice with the same propertyll 
Compute the dual lattice 

A = jxVAfo; rfieV) (15) 



with • A^ = <5^. The basis vectors \%, have the obvious property that A*^ • 

© Now try to build linear combinations 

tA := Y. mPX h 

P 

such that 

t A ■ g (^W) = l mod 2 . (16) 

It is obvious that one just has to scan all possible combinations with mP G {0, 1} 
since an even will always lead to an even number on the right-hand side of 
( fTBl) . tA is then unique up to U(l) generators under which all ip^ 1 ' are neutral. 

© Given a generator tA, one has to check whether the model contains fields <f>w with 
charges of the type 'even over odd'. If this is the case, switching on the fields 
yields configurations with a Z2 symmetry under which the ip® fields are odd. 

Only the last two steps have to be slightly modified in order to obtain an arbitrary Zjv 
symmetry. In section 15.21 we will apply these methods in order to identify phenomeno- 
logically attractive string vacua. 



2 This lattice can, for instance, be obtained with the Mathematica command 'LatticeReduce' 
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3 Simplifying multiple Zjv symmetries 



In this section we will consider a finite Abelian group G = x . . . x Z^ with K fields 
z = 1, . . . , K, transforming under G. This setup may or may not be a result of the 
diagonalization procedure described in section [2][f| Our aim is to eliminate redundancies, 
i.e. make the discrete symmetry 'as simple as possible'. 

First, consider a toy example which consists of a Z^-symmetry with one field ip with 
charge 4. This is equivalent to a Z^-symmetry where ip has charge 2. In the case of 
one Z^r-symmetry with fields ip® and charges q(ip®), it is easy to see if there is some 
redundancy in our description of the symmetry. If the greatest common divisor (GCD) 
of the q(ip®) and N is greater than one, we can divide q{ip®) and N by their GCD. 
That is precisely what happened in our above toy example. This idea of reducing the 
group order carries over to the general case. In addition, the method we are going to 
develop will bring G into a canonical form. Every finite Abelian group can be written 
SIS db direct product of the form 

Canonical = Z dl X . . . X Z djv , where di divides d i+1 . (17) 

The di are uniquely determined by the group and are called invariant factors. The 
diagonalization process described in section [2] always leads to this form (cf. equation ((7j)). 



3.1 The general case 

Let us look at the charge matrix (Q^ij = qj(ip®). This is a K x N matrix with integer 
entries. Again, the idea is to calculate the Smith normal form of Q^. However, to get a 
meaningful result, we need to bring G into the form Z^ x . . . x Z^ where d is the least 
common multiple (LCM) of the di. This is because a Z^ discrete symmetry allows us to 
perform discrete rotations of the generators. Enlarging the symmetry implies a rescaling 
of the charge matrix, (Q'^p)ij = j- qj(ip®). The fields transform now according to 



exp ^27ri ^{Q'^a^j 



?/> W -> exp ( 27ri 22(Q^)ij % J where a 6 = -j- with < l t < d-1 . (U 



(QL) can be brought into Smith normal form S by unimodular transformations E G 
GL(AT,Z) and F G GL(iV, Z), 

EQ'^F = S = diag'(si, . . . , Sfc) where k = mm(K,N). (19) 

If rankQ^ < k, some Sj might vanish. The corresponding rows do not have to be 
considered further. This yields the transformation behavior 



^ - exp 2th E V Sim F~l a n j,® 

= exp^Tri YEj 3j F7^) (20) 



3 That is, we do not require di\di+\ here. 
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Just like in equation (jHJ), we are allowed to define 



y F- 1 



Hence, we get 



(i) 



exp hvri^^ 1 



d/sj 



(0 



(21) 



This tells us that we have rewritten our symmetry group G as Z^ x . . . x Z<j' where d\ 
is the numerator of the (reduced) fraction j- and d' i+1 divides d\. If a d\ is equal to one, 
this factor can be omitted. The new charges of the fields, q"(tp^), are encoded in the 
matrix E" 1 , 



E~ 



( q"{ip {1) ) \ 
q »{^)) 



\ 



(22) 



/ 



These charges are equivalent to the the discrete charges s[ q'{^^) (with denoting the 
denominator of the reduced fraction d/si), which one immediately reads off. An obvious 
consequence of our discussion is that after simplification there are at most as many 
Z<^ factors as fields. Note that the volumes of the fundamental region of G"s lattice 
{Vg = YliLi di) and the lattice spanned by the ip® charges (V^ = detQ^) provide a 
necessary but not sufficient criterion for redundancies: if both possess a GCD, the size 
of the symmetry may be reduced by this GCD. 



3.2 An alternative derivation 

Before discussing an example, let us present an alternative point of view. A coupling 
(■0W) m . . . (■0C f O) n -K is allowed by the discrete symmetry G = Z^ x . . . x Z^ only if 
there is a vector n e Z x such that 

Q^n = diag(c/i, . . . , d N ) m (23) 
with some m G Z N and Ql = (g(^ (1) ), • • ■ ,q(^ (K) ))- Equation (J23D can be rewritten as 

diag {i--£)- Q > = dm > (24) 

where d denotes the LCM of the d{, as before. Now we diagonalize the matrix on the 
left hand side of the equation, 

diag -^VqJ = (F-Y-S-iE-y (25) 



d\ d 



N 



with the N x K matrix S = diag'(si, . . . Sk), the unimodular matrices E and F, and 
k = mm(K, N), as before. Let v denote the rank of Q^, i.e. S = diag'(si, . . . s u , 0, . . . ). 
Now (l23l) can be recast as 



diag'(|,...,^,0,... ) (E-Yn = m>, (26) 
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where w! = F T m, and again and d\ denote the numerators and denominators of the 
reduced fractions Si/d, respectively. The rank of the matrix on the left-hand side of this 
equation is v. We can therefore truncate the equation, 

diag(|,...,^) {E'Yn = m u , (27) 

where {E' 1 ) denotes the left v columns of E~ x (such that (E~ l ) T is the upper v x K 
part of (E~ l ) T ), and m v G 7L V . This equation is equivalent to 

{E-^n = diag' ( % . . . M m v , (28) 

Comparing this equation with (1231 reveals that the rows of E~ x contain the charges of 
the ip^ and the d\ determine the canonical symmetry. 



3.3 Example 

Let us continue the example of section [2731 After breaking both U(l)s we are left with 
a 7L>2 x 7L§ and a charge assignment given in table [2j We have to extend our symmetry 



(a) Original (b) 'Blown up' (c) Minimal 

charges. charges. charges. 

1 1 3 1 ipW 1 

^ 1 3 V (2) 3 3 ^ (2) 3 



Table 2: The example from section 12.31 continued. The original charges (a) are blown 
up to the charges of an extended Z 6 x Z 6 symmetry, which can be reduced to a Z 6 
symmetry by discrete rotations. 

to 7jq x The charge matrix and the Smith normal form are 

The d' t can be inferred from the diagonal matrix S: they are given by 6 times the inverses 
of the diagonal entries, i.e. we have d[ = 6 and d' 2 = 1. The numerators are hence 6 and 
1, such that we are left with a Zj x Zi = symmetry. The charges are given by the 
rows of E^ 1 (modulo 6); since one factor is trivial we obtain that ip^ has charge 1 and 
ip( 2 ^ has charge 3. 
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e 2 



^(2) 




^,(2) 



(a) Extended lat- 
tice. 



(b) Simplified lat- 
tice. 



Figure 2: (a) Extended and (b) simplified charge lattices. 



3.4 Visualization 

The visualization works as before. We extend the lattice to TL§ x 7jq (figure [2] (a)). 
Our simplification process amounts to identifying a direction on which both nodes lie 
(figure [5] (b)). It is easy to see that ift^ 1 ' sits at 1/6 of the length of the one- dimensional 
7L§ lattice, corresponding to Iiq charge 1, while ip^ sits at 1/2 = 3/6 of the length, which 
leads to Z 6 charge 3. 



4 Automatization 

We provide a Mathematica package which automatically identifies the remnant symme- 
tries, as discussed in section |2} and brings them to the canonical form, as described in 
section [3], on our web site [11]. This package has been used in our applications, which 
will be discussed in what follows. 



5 Applications 

5.1 GUT model building 

Let us now apply the above methods to model building. We focus on a specific GUT 
model, which has been discussed in [12]. There, an anomaly-free 7jq symmetry was 
found that may allow us to suppress proton decay in SO(10) GUTs. The field content is 
given in table [31 There are three generations ip m of standard model matter; H and H' 

</> m H H' tfj H A S 

SO(10) 16 10 10 16 16 45 54 

Z 6 ' 1 -2 2 -2 2 

Table 3: Field content of an SO(10) model. 



10 



contain the standard model Higgses. The other fields are used to break SO(10) down to 
the standard model gauge group Gsm- Since SO(10) has rank 5 there is a U(l) factor, 
called U(l) x , in addition to the standard model gauge group. This U(l) x gets broken by 
giving VEVs to the SM singlet fields contained in ipu and ip H , which have U(l) x charges 
±5. Hence, we have a situation where a U(l) x x 1, e symmetry gets broken by a field 
with charges (±5, ±2). We can apply the routine described in sectional This can be 
accomplished by extending the to a U(l), introducing a dummy field with charges 
(0,6) and assigning a VEV to this field. 

We find that the U(l) x x Z 6 gets broken to a Z 30 with charges given in tablelU Let us 

Q U D L E H u H d 
Z 30 3 3 11 11 3 24 16 
Z 5 xZ 6 (1,1) (1,1) (2,5) (2,5) (1,1) (3,4) (2,0) 

Table 4: The Z 30 -Z 5 xZ 6 charges of the MSSM field content. 

remark that the Z5 subgroup of this Z30 is redundant in the following sense: whenever 
a coupling is gauge invariant under Gsm ; the coupling is also invariant under Z5. This 
is because SU(5) has a non-trivial center Z 5 (cf. the analogous discussion in [13]). Note 
also that Gsm invariance is the same as SU(5) invariance since the Cartan generators 
are equivalent. Our results imply that, if we assume that all states except the ones 
listed in table H] attain masses, we are left with an anomalous field content. Specifically, 
the three generations of SM matter are anomaly free, but the 7jq charges of the Higgs 
fields exhibit an anomaly (for a discussion of discrete anomaly constraints see [14-16]). 
Using discrete anomaly matching [13] we can hence infer that the above assumption is 
inconsistent: either further light states have to be present, or the 7*6 cannot be exact, i.e. 
we have to introduce further fields that attain VEVs. Of course approximate symmetry 
might be sufficient for suppressing the dangerous dimension five operators, and might 
well be correlated with flavor hierarchies (cf. the discussion in [17]). On the other hand, 
our findings show that the TL§ symmetry introduced in the SO(10) GUT in [12] cannot 
give rise to proton hexality [10]. 

5.2 String model building 

As mentioned in the introduction, an exact matter parity can also be obtained in string 
theory by breaking a U(1)#-l symmetry by two units [4,18,19]. This has led to a 
couple of vacuum configurations with an exact matter parity. Yet it turns out that, 
with this strategy, one is not always successful: within the so-called mini-landscape [18] 
of heterotic orbifolds with exact MSSM spectra, vacuum configurations with an exact 
matter parity could only be identified in a small fraction of the models; that is, an 
appropriate U(l) b-l symmetry could only be identified in 15 out of 218 possible models. 
The obstacles encountered in this U(l)s_i,-based approach are perhaps best illustrated 
in a concrete example. The model discussed in [20] (which later became absorbed in 
the mini-landscape) does have a U(1)b-l symmetry [19,21], yet in the 4D zero mode 
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spectrum there is no field with an even B — L charge (nor with a fractional charge of 
the type 'even over odd', which is also sufficient to provide us with matter parity [4])0 
On the other hand, it is also clear that one does not really need a U(l) b-l symmetry 
with the standard charges for the MSSM matter fields. Any U(l) symmetry under 
which the matter states have odd charges, and for which there exist SM singlets with 
even charges, could do the job. Using the methods discussed in section l2~5l we were able 
to identify a collection of Gsm invariant fields that break the U(l) factors down to 
matter parity in the model presented in [20,21]. It is given by 

{0 W } = {Si, S 2 , S 3 , S 5 , S 7 , Sg, S U , S M , S W , S 18 , S 19 , S 2 0, S 2 2, S 2 3, «24, S34, S39, S40, 

S41, S48, S53, S54, S57, S58, S59, S60, Sqi, Sq2, $65, s 66> fl-i, /l-4j ^1-14 } (30) 

in the notation of [21]. Our algorithm gives us a Zio symmetry, however, as discussed in 
section EIH the Z 5 , which is just the non-trivial center of SU(5), is redundant. A quick 
scan indicates that in many (if not in all) mini-landscape models vacua with matter 
parity can be obtained. A detailed analysis of these issues and of the phenomenological 
properties of such vacua will be carried out elsewhere. 

6 Discussion 

We have described a simple method to determine symmetry breaking patterns which 
arise when U(l) gauge theories get broken to discrete subgroups. This method has a 
very simple geometrical interpretation: the fields acquiring VEVs define a charge lattice. 
Couplings (ip^) ni (-f^ 2 ))™ 2 . . . are only allowed by the remnant discrete symmetries if the 
sum of the charge vectors, niqfy^) +n2q{ip^) + ■ ■ ■ , lies on a node of the charge lattice. 

Unimodular transformations allow us to identify the remnant discrete symmetries, 
and to make them as simple as possible, i.e. to determine the true (or canonical) sym- 
metries in a unique way. 

We have applied our methods to model building. In the context of GUTs, we have 
identified an obstacle to completely forbidding dimension five proton decay operators in 
certain SO(10) GUTs. In string model building our methods allow us to identify novel 
vacuum configurations with an exact R parity. 
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